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ABSTRACT 


The  annual  cycle  of  the  North  Pacific  wind  stress,  Ekman  pumping  and  Sverdrup 
transport  is  investigated  by  means  of  empirical  orthogonal  function  (EOF)  analysis 
techniques.  Fifty-two  years  of  National  Centers  for  Environmental  Prediction  (NCEP, 
formerly  "NMC")  and  the  National  Center  for  Atmospheric  Research  (NCAR)  Reanalysis 
daily  averaged  surface  wind  components  covering  the  extratropical  North  Pacific  are 
used  to  calculate  daily  averaged  wind  stress  components.  These  wind  stress  components 
are  averaged  to  624  monthly  mean  fields  from  which  monthly  mean  Ekman  pumping  and 
Sverdrup  transport  fields  are  derived.  Each  data  field  has  the  long-term  annual  mean  and 
linear  trend  removed  at  each  grid  point  before  EOF  analysis.  The  first  three  modes  are 
considered  to  be  physically  significant  based  on  the  plot  of  the  log  eigenvalue  against  its 
root.  The  first  principal  components  of  each  field  are  highly  coherent,  have  most  of  their 
spectral  energy  at  the  annual  cycle,  and  are  nearly  in  phase.  The  second  mode  of  these 
fields  show  significant  spectral  energy  at  the  semiannual  cycle.  The  third  mode  has  more 
complex  spatial  variability  than  the  first  and  second  mode.  The  annual  cycle  is  best 
described  by  the  first  EOF  mode  which  is  composed  of  two  seasons  (winter  and  summer) 
that  last  for  five  months  each  separated  by  two  rapid  (approximately  one  month) 
transition  periods  (spring  and  fall).  The  annual  cycle  is  complex  and  certain  geographical 
regions  are  seasonal  "hot  spots"  that  add  to  the  complexity.  Relationships  are  derived 
between  the  various  modes  to  help  explain  seasonal  movement  of  the  Ekman  pumping 
zero  isotach  and  seasonal  variations  in  Sverdrup  transport.  The  annual  variability  can 
reasonably  be  described  using  only  the  first  three  EOF  modes. 
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I.  INTRODUCTION 


The  surface  momentum  flux  from  the  atmosphere,  or  wind  stress,  and  the  curl  of 
the  wind  stress  acting  on  the  sea  surface  are  fundamental  forcing  agents  for  dynamic 
ocean  processes  (Bakun  and  Nelson,  1991).  Basin-wide  winds,  rather  than  local  forcing, 
have  been  shown  to  be  responsible  for  interannual  fluctuations  in  the  California  Current 
System  (CCS)  (Chelton  and  Davis,  1982).  Bryden  et  al.  (1991)  suggest  that  the  surface 
circulation  in  the  North  Pacific  controls  the  net  meridional  heat  transport.  Large-scale 
wind  stress  variations  and  the  resulting  Ekman  pumping  have  major  effects  on  the 
biological  systems  and  marine  fisheries  at  El  Nino  and  decadal  time  scales  (Brodeur  and 
Ware,  1992;  Beamish  and  Bouillon,  1993;  Parrish  et  al,  2000).  Ocean-atmosphere 
coupled  models  suggest  that  ocean  temperature  variability  is  due  to  coupling  between  the 
subtropical  gyre  and  the  wind  stress  field  associated  with  the  Aleutian  Eow  (Eatif  and 
Barnett,  1994,  1996). 

There  are  no  direct  long-term  measurements  of  ocean  circulation  with  regular 
spatial  and  temporal  sampling.  Circulation  fields  can  be  estimated  from  modeled  and 
observed  wind  data  by  using  fluid  dynamical  theory  to  derive  lateral  and  vertical  currents 
(Ekman,  1905;  Sverdrup,  1947;  Stommel,  1948;  Munk,  1950).  In  recent  years,  ocean 
circulation  estimates  based  on  wind  stress  data  from  either  observations  (Wyrtki  and 
Myers,  1976,  Hellerman  and  Rosenstein,  1989,  Chelton  et  al.,  1990,  and  Xie  and  Hsieh, 
1995)  or  model  output  have  been  developed  for  individual  ocean  areas  and  for  the  world 
ocean  (Trenberth  et  al.,  1990). 

The  goal  of  this  thesis  is  to  describe  the  annual  cycle  of  the  extratropical  North 
Pacific  wind  stress  and  the  derived  circulation  fields  (Ekman  pumping  and  Sverdrup 
transport)  using  Empirical  Orthogonal  Eunction  (EOE)  analysis.  Annual  means  of  the 
output  of  the  EOE  analysis  and  statistical  correlations  between  the  modes  of  the  various 
fields  will  be  described  in  an  attempt  to  identify  their  seasonal  dominance  and  spatial 
patterns.  Table  1  provides  a  list  of  symbols  used  and  follows  standard  oceanographic  and 
statistical  practice. 
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Variable 

Deseription 

Units 

u  (u) 

Daily  mean  eastward  (northward)  wind  stress 
eomponent. 

Pa 

Monthly  mean  eastward  (northward)  Ekman  transport 

2  -1 

eomponent. 

m  s 

UTT 

Monthly  mean  eastward  (northward)  wind  stress 

Pa 

eomponent. 

W 

Ekman  pumping  vertieal  veloeity. 

em  day'' 

(r) 

Monthly  mean  Sverdrup  transport. 

10'’  m^  s'' 

Pa 

Density  of  air 

1 .22  kg  m'^ 

<P 

Geographie  latitude 

Degrees 

X 

Geographie  longitude 

Degrees 

X 

Input  data  matrix  for  a  single  month 

X 

Entire  time  series  of  the  input  data  matrix 

X  =  T,w,if/ 

X 

Eong-term  annual  mean  of  input  data  matrix 
(detrended) 

N 

Number  of  input  matrix  grid  points. 

M 

Number  of  months  in  the  data  time  series. 

624 

n 

Grid  point  index 

n=(l,2,...N) 

m 

Time  index 

m=(l,2,...M) 

PCt 

k^*’  prineipal  eomponent  of  X 

k=(l,2,...N) 

E" 

k^*’  eigenveetor  of  X 

k=(l,2,...N) 

k'*’  eigenvalue  of  X 

k=(l,2,...N) 

Mf 

k*'’  mode  of  X  refers  to(  PC^ ,  )  pairings 

H 

Varianee/eovarianee  matrix  of  X 

m 

Subseript  denoting  the  arithmetie  mean  for  eaeh 

1  ^ 

rn  =  —  /  m.: 

ealendar  month  (m)  over  the  entire  time  series  (Y  = 

(624  months  x  year/12  months)  =  52  Years 

1=1 

=  (1,2,. ..12) 

Table  1 .  Variables  Deseription  and  Units  Used  in  Text. 
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II.  DATA  AND  METHODS 


A.  DYNAMICAL  CALCULATIONS 

A  52-year  subset  (1948-1999)  of  daily  averaged  surfaee  wind  eomponents  from 
the  National  Centers  for  Environmental  Predietion  (NCEP)  and  the  National  Center  for 
Atmospherie  Researeh  (NCAR)  Reanalysis  projeet  (Kalnay  et  ah,  1996)  eovering  the 
North  Paeifie  was  used  to  ealeulate  daily  averaged  wind  stress.  The  geographieal  limits 
used  are  15°N  to  70°N  and  100°E  to  80°W.  The  data  were  gridded  at  2.5°  intervals  for 
both  latitude  and  longitude. 

Daily  averaged  wind  stress  eomponents  were  ealeulated  using  the  formula: 


w. 


995 


u, 


995  ’ 


w. 


995 


995 


(1) 


where  (r^)  denotes  the  eastward  (northward)  daily  averaged  eomponent  of  wind 
stress,  is  the  density  of  air  (a  eonstant  1 .22  kg  m'  ),  is  a  non-linear  drag  eoeffieient 
based  on  Large  and  Pond  (1981)  and  modified  for  low  wind  speeds  as  in  Trenberth  et  al. 
(1990),  IW995I  is  the  wind  speed,  and  U995  (V995)  is  the  eastward  (northward)  eomponent 
of  the  wind  veloeity  (m  see'')  at  sea  level,  respeetively. 

The  resultant  long-term  monthly  mean  eastward  (T^)  and  northward  (T^)  wind 
stress  eomponents  at  eaeh  station  (N  =  1679  grid  points)  were  eomputed  as  the  arithmetie 
means  of  the  daily  eastward  and  northward  eomponents  for  eaeh  month  over  the  entire 
time  series  (M  =  624  months).  These  values  were  used  to  derive  the  monthly  mean 
eireulation  fields  as  follows. 

The  monthly  averaged  eastward  (M^)  and  northward  (M^,)  eomponents  of  the 
2  1 

Ekman  transport  (m  s'  )  were  ealeulated  using  the  formula: 
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(m.,  m,,) 


T, 


T 


(2) 


V  Pswf  Pswf ^ 


where  is  the  density  of  seawater  (a  eonstant,  1025  kg  m'  )  and  f  is  the  Coriolis 
parameter.  Monthly  averaged  values  for  Ekman  pumping,  or  open-oeean  upwelling  (em 
day"')  were  ealeulated  as  the  divergenee  of  the  surfaee  Ekman  transport  using  the 
formula: 


^iij)  = 


/'  M  -M  ^ 


+ 


2A^ 


(3) 


where  (j)  and  X  denote  geographie  latitude  and  longitude,  respeetively,  AA,  =  2.5°  and  Acj)  = 
(AA,  eos  (j)).  Subseripts  i  and  j  represent  eolumn  and  row  position  in  the  data  matrix 
(Eigure  1). 


ij+l 

• 

•  • 

T 

A^ 

ij  i+ki 

• 

•  ^ —  A/i  — ^  • 

• 

G./-1 

•  • 

Eigure  1 .  Grid  Spaeing  Nomenelature  Eigure.  The  Variables  /  And  J  Represent  Column 
and  Row  Position,  Respeetively,  in  the  Data  Matrix.  <!>  and  X  Denote  Geographie 
Eatitude  and  Eongitude  Respeetively. 
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Centered  differenees  were  used  throughout  the  interior  of  the  data  matrix. 
Column  and  row  end  values  were  used  at  the  boundaries  of  the  data  matrix.  Most  of  the 
eolumn  and  row  end  values  lay  over  landmasses  and  any  artifieial  feature  introdueed  at 
the  boundaries  did  not  greatly  affeet  oeeanie  eireulation  ealeulations. 

Sverdrup  transports  were  eomputed  from  the  meridional  transport,  .  This  was 

derived  from  the  Sverdrup  transport  streamfunetion  relationship  by 

0 

integrating  westward  from  the  eastern  oeean  boundary  to  the  western  oeean  boundary 
between  15°  and  57.5°N  (N  =  1314  grid  points).  The  Aleutian  Islands  were  ignored. 

6  3  1 

Monthly  averaged  values  for  Sverdrup  transport  Sverdrup  units  =10  m  s'  )  were 
ealeulated  by: 


AA 


(m  +M  ,  , 


(4) 


where  is  the  Sverdrup  meridional  transport  at  loeation  and  is  initialized  to 

zero  at  the  eastern  boundary. 

B,  ANNUAL  MEAN  FIELDS 

Annual  mean  fields  were  eomputed  as  the  arithmetie  mean  over  the  entire  time 
series  and  are  shown  in  Figure  2.  The  mean  wind  stress  field  (Figure  2a)  is  dominated  by 
a  eyelonie  eireulation  eentered  about  52°N,  170°E  assoeiated  with  the  Aleutian  low  and 
antieyolonie  eireulation  assoeiated  with  the  subtropieal  North  Paeifie  high  eentered  about 
30°N,  140°W.  There  are  strong  westward  wind  stresses  in  the  subtropies  (15°-25°N)  and 
strong  eastward  wind  stresses  in  the  mid-latitudes  (35°-50°N,  140°E  -  180°).  The 
eastward  wind  stress  drives  a  southward  Ekman  transport  earrying  eold  water 
equatorward.  The  persistent  westward  trade  winds  in  the  subtropies  drive  a  northward 
Ekman  transport  of  warm  water  poleward.  Due  to  eonvergenee  (divergenee)  of  the 
surfaee  Ekman  layer  within  the  westerly  (easterly)  wind  stress  regions,  the  subpolar 
region  has  resultant  open  oeean  upwelling  and  the  subtropieal  region  has  resultant  open 
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ocean  downwelling  (Figure  2b).  Strong  equatorward  alongshore  and  divergent 
(upwelling  favorable)  wind  stress  is  found  along  California  and  the  Baja  peninsula.  The 
equatorward  wind  stress  along  the  eastern  boundary  ereates  a  region  of  intense  coastal 
upwelling.  Through  conservation  of  potential  vortieity,  antieyelonie  cireulation  in  the 
subtropical  gyre  and  cyclonic  circulation  in  the  subpolar  gyre  is  established  (Figure  2c). 


Figure  2.  Mean  Fields  Computed  for  the  Period  1948-1999.  (a)  Wind  Stress.  Scaling 
Vector  is  0.25  Pa.  (b)  Ekman  Pumping.  The  Contour  Interval  is  5  cm  day'\  Zero  Isotach 
is  Bold.  Upwelling  Regions  are  Shaded,  (c)  Sverdrup  Transport.  The  Streamline 
Interval  is  10  Sv  (1  Sv  =  10  m  s'  ).  Arrowheads  Indicate  the  Direetion  of  the 
Transport.  Zero  is  Bold.  Positive  Values  are  Dashed  Lines  and  Negative  Values  are 

Solid  Lines. 
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Figure  3.  Standard  Deviation  of  the  Mean  Fields  Computed  for  the  Period  1948-1999.  (a) 
Wind  Stress.  Sealing  Veetor  is  .07  Pa.  (b)  Ekman  Pumping.  Contour  Interval  is  5  em 
day'  .  (e)  Sverdrup  Transport.  The  Streamline  Interval  5  Sv  (1  Sv  =  10  m  s'  ). 
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C.  EOF  ANALYSIS 

Atmospheric  forcing  and  the  resultant  oceanographic  response  are  variable  on  a 
range  of  temporal  and  spatial  scales  (Hare,  1996).  These  overlapping  processes  produce 
a  complicated  picture  that  often  defies  simple  explanation.  However,  with  the  use  of 
"eigentechniques",  such  as  EOF  analysis,  one  can  link  both  the  temporal  and  spatial 
patterns  of  long-term  records  covering  large  spatial  areas  (Kutzbach,  1967;  Barnett,  1977; 
Hardy  and  Walton,  1978;  Legler,  1983;  Hare,  1996;  Mizoguchi  et  al.,  1999). 

The  underlying  mathematical  principles  of  EOF  analysis  for  scalar  values  (Ekman 
pumping  and  Sverdrup  transport)  and  vector  quantities  (wind  stress)  are  similar.  The 
goal  of  such  an  analysis  is  to  reduce  the  number  of  variables  in  the  original  data  set  while 
maintaining  the  total  variance  of  the  original  data  set  (Hare,  1996). 

The  EOF  analysis  used  in  this  study  parallels  the  notation  and  methodology  of 
Eegler  (1983).  A  data  field  was  defined  over  N  grid  points  on  M  time  steps.  Each  data 
field  had  the  long-term  annual  mean  (X)  and  linear  trend  removed  at  each  grid  point. 
The  input  matrix  X  was  of  dimension  N x  M,  with  n=(l,2,  ...N)  indexing  grid  points  and 
m=(l,2,...M)  indexing  time,  and  contained  the  whole  data  set.  The  value  of  the  field  at 
any  point  in  time  and  space  is  denoted  as  x„_m-  Wind  stress  components  at  each  station 
were  combined  into  a  vector  quantity  of  the  complex  form  x„  „  =  T^  +  i  T„  .  The 

input  data  matrix  is  written  in  matrix  notation  as 


Xji  Xj2 
^21 


x  = 


^IM 


(5) 


^N\ 


^NM 


The  variance/covariance  matrix  (//)  is  expressed  as 


8 


(6) 


where  represents  the  transpose,  or  complex  conjugate  transpose  in  the  case  of  the 
wind  vectors,  of  X.  //  is  a  Hermitian  matrix,  which  is  symmetric  (NxN)  with  real¬ 
valued  diagonal  elements. 

The  EOF  analysis  has  three  results:  principal  components  (PCf),  eigenvectors 
{Ef)  and  eigenvalues  (/if),  where  the  subscript  k=  (1,2,...N)  and  the  superscript  X 
designates  the  data  field  (wind  stress  (r),  Ekman  pumping,  (w),  and  Sverdrup  transport 
{y/).  The  eigenvalue  decomposition  of  [//]  produces  a  diagonal  matrix  [/if]  of 

eigenvalues  and  a  full  matrix  [Ef^]  whose  columns  are  the  corresponding  eigenvectors  so 
that 


(7) 


Each  eigenvalue,  /if,  has  an  associated  spatial  eigenvector,  Ef .  The  eigenvectors 
have  a  spatial  orthogonal  relationship,  that  is,  each  eigenvector  is  orthogonal  to  all  those 
preceding  it.  are  in  the  complex  form  (Re  (Ef )+  ilm^Ef^). 

Ef  and  Ef  fields  are  displayed  as  contour  and  streamline  maps,  respectively.  El 
fields  are  displayed  in  vector  form.  El  vectors  are  computed  using  the  formula 

(£;)=  (8) 


The  PC  for  each  eigenvector  k  at  time  m  is  obtained  from  the  equation 
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(9) 


PC  = 

^  ^kfh 


(£/)'  DK.] 


where  X-  represents  a  column  vector  from  the  monthly  mean  of  the  input  data.  The  time 
series  of  the  coefficients  for  eigenvector  Ef  is  known  as  the  k‘’'  principal  component 
{PCf  ).  In  this  study  a  seasonal  time  series  of  PCs  are  computed  as  the  arithmetic  mean 
(m  )  of  all  the  PCs  for  each  month  (m)  over  the  entire  time  series  m=(l,2,  ...M). 

The  spatial  positioning  of  the  eigenvector  patterns  does  not  change  over  time. 
PC'^  modulate  the  magnitude  of  the  eigenvector  fields  over  time.  PC]^  modulate  both 
vector  magnitude  and  rotation  angle  over  time.  PCI  is  of  the  complex  form 
Re  (PCl-  )+  i  Im  (PC^^.  Magnitude  and  rotation  angle  information  are  calculated  as: 


magnitude  [PCI^ )  =  ^Re(RC;^)"  +Im(RC;^)" 


and 


rotation 


angle  [pCl-)=i&n  ^ 


(10) 


(11) 


El  vector  fields  are  modulated  by  both  the  magnitude  {PCl^ )  and 
rotation  angle  [pCl-)  time  series.  A  positive  (negative)  rotation  angle  [pCl-]  rotates 
El  vectors  in  a  cyclonic  (anticyclonic)  manner  around  the  grid  point  at  the  base  of  each 
vector. 
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The  PC^,E^  pairings  are  referred  to  as  modes  and  are  designated  Mf .  The 
modes  are  ranked  by  their  overall  variability  with  respeet  to  the  total  dataset.  The 
eontribution  of  eaeh  mode  to  the  total  varianee  of  the  dataset  is  ealeulated  by  dividing  the 
eorresponding  Zf  by  the  traee  of  the  varianee/eovarianee  matrix  where 

N 

TRACE  (12) 

k=\ 


The  original  data  at  any  time  m  can  be  recreated  by  the  complete  set  of 
orthonormal  eigenvectors  through  the  relationship 


N 

V  i=\ 


(13) 


It  is  important  to  remember  in  the  following  diseussion  that  the  output  of  this  EOF 
analysis  represents  anomalies  to  the  mean  field. 
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III.  RESULTS 


A.  STATISTICAL  ANALYSIS 

Two  methods  were  used  to  determine  the  signifieanee  of  the  EOF  modes  that 
were  based  on  the  relative  magnitudes  of  the  eigenvalues;  the  Seree  test  and  the  Log 
Eigenvalue  (LEV)  graph.  The  term  seree  derives  from  the  resemblanee  to  the  rubble  that 
forms  at  the  foot  of  a  mountain.  Cattell  (1966)  hypothesized  the  seree  represented 
unwanted  noise  and  that  only  the  EOF/PCs  prior  to  the  seree  should  be  retained  for 
further  use.  Identifieation  of  the  exaet  root  where  the  scree  begins  to  form  is  rather 
subjective,  however.  A  visual  inspection  of  a  plot  of  the  eigenvalue  versus  its  root  (Scree 
test)  (Figure  4)  showed  the  typical  eigenvalue  plot  that  had  a  steep  slope  over  the  first 
few  roots  and  then  a  gradual  trailing  off  of  the  rest  of  the  roots  (the  scree).  The  first  two 
eigenvalues  define  the  slope  and  the  third  is  at  the  top  of  the  scree.  In  other  words, 

appear  to  be  above  the  scree  and  therefore  may  be  significant. 


Figure  4.  Scree  Test  Plot  of  /I  [  (Solid  Line),  (Solid  Bold  Line)  and  X  ^ /lO^^ 

(Broken  Bold  Line)  for  k=(l,2,. .  .10)  versus  Eigenvalue  Root  Number. 

A  plot  of  the  LEV  versus  its  root  (Craddock  and  Flood,  1969)  is  shown  in  Figure 
5.  Experimental  results  indicate  that  an  EOF  of  random  data  results  in  a  straight-line 
LEV  graph  (Farmer,  1971).  To  determine  the  number  of  modes  to  retain,  a  straight  line 
is  drawn  through  the  higher  number  roots,  and  those  lying  above  the  line  are  retained.  To 


13 


emphasize  the  break-point  of  the  lower  number  roots  from  the  straight  line  (gray  line), 
only  the  first  400  modes  are  shown  here.  were  signifieantly  above  the  noise  level 

and,  therefore,  may  have  some  physieal  meaning. 


Figure  5.  Log  Eigenvalue  (LEV)  Graph  for  the  Data  Eields  as  Labeled.  The  Eirst  Three 

Eigenvalues  are  Designated  by  Triangles. 

This  study  foeused  on  M^’^f  .  are  presented  in  Appendix  D.  The  pereent 

of  varianee  and  eumulative  varianee  for  the  first  three  eigenvalues  are  given  in  Table  2. 
Only  a  few  EOE  modes  were  needed  to  explain  a  large  part  of  the  varianee  for  eaeh  series 
of  fields.  V/j^2  3’  M'23’  aeeounted  for  57.9%,  38.5%  and  74.6%  of  the  total 

varianee,  respeetively. 
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VARIANCE 


Mode 

r 

w 

¥ 

Percent 

Cumulative 

Percent 

Cumulative 

Percent 

Cumulative 

1 

34.5 

34.5 

23.5 

23.5 

45.6 

45.6 

2 

14.9 

49.4 

9.9 

33.4 

19.7 

65.3 

3 

8.5 

57.9 

5.1 

38.5 

9.3 

74.6 

Table  2.  Percentage  of  the  Total  Variance  and  the  Cumulative  Variance  Explained  by  the 

2  2 
Eirst  Three  EOE  Modes  for  Wind  Stress  (  r  =  Pa  ),  Ekman  Pumping  {w  =  cm 

2  2 

day'  )  and  Sverdrup  Transport  =  Sv  ). 

Spectral  analysis  of  the  entire  time  series  of  (Eigures  7a,  b,  c)  showed  that 

most  of  spectral  energy  was  found  at  the  annual  (0.083/month)  cycle.  PC^’'^’’^ 

also  exhibited  a  small  amount  of  energy  at  the  semiannual  (0.167/month)  cycle. 
PCl’'^  spectra  had  more  energy  at  the  semiannual  cycle  than  at  the  annual  cycle;  however 
PC2  spectra  had  more  energy  at  the  annual  cycle  than  at  the  semiannual  cycle.  Power 
spectra  for  PCl’'^  was  generally  in  the  noise  level  of  the  spectra,  although 
ECj"  spectra  showed  energy  above  the  noise  level  at  the  semiannual  cycle.  Coherence 
function  estimates  for  the  entire  series  of  PC^’'^''^  (Eigure  7)  showed  high  coherence  at  the 
annual  frequency. 
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month 


ire  6.  Power  Speetra  of  PC^  (Thin  Solid  Line),  PC2  (Broken  Line),  and  PCf  (Bold 
Solid  line)  for  (a)  Re  {pCl  \  (b)  Re  (rC^)  and  (e)  Re  {PCl )  for  the  Entire  Time  Series. 
Eaeh  time  Series  has  been  Weighted  by  its  Varianee. 


Figure  7.  Coherence  Function  Estimate  for  the  Entire  Monthly  Averaged  Time  Series 
(1948-1999)  of  Re(ECi")  and  PC'^  (Solid  Eine),  Re(EC;)  and  PC'[  (Broken  Eine), 

PC^  And  PC^  (Broken  Dotted  Eine). 


B,  EOF  MODE  1 

Ml  (Figures  8a,  b  and  9a),  (Figures  8c  and  9b)  and  Ml'  (Figures  8d  and  9c) 
account  for  34.5%,  23.5%  and  45.6%  of  the  total  variance,  respectively  (Table  2).  The 
first  principal  components  for  wind  stress,  Ekman  pumping  and  Sverdrup  transport, 
designated  PC^^ ,  were  uniquely  similar  in  that  all  three  have  a  coherent  annual  cycle. 

PC^-  had  an  annual  cycle  with  a  maximum  in  January  and  a  minimum  in  July,  August 
and  September  for  Re(E’Qc),  ECj^and  PC^-,  respectively.  E’Ci'k  had  a  larger 
magnitude  during  the  winter  maximum  (Re(E’Cj),,)  =  1.375,  PCl"^  =  0.059  and  = 
0.052,  where  m  =  January)  as  opposed  to  the  summer  minimum  (Re(ECj''^)  =  -1.135, 
PCl'm  =  -0.051  and  PCl^„  =  -0.035,  where  m  =  July,  August  and  September, 
respectively).  The  range  of  one  standard  deviation  was  also  greater  during  the  winter 
months  as  opposed  to  the  summer  months.  In  other  words,  the  wind  stress  was  stronger, 
pumping  more  vigorous,  and  the  transport  more  robust  during  the  height  of  winter  than 
during  the  peak  of  summer.  There  was  also  greater  directional  and  strength  variability 
among  all  these  fields  during  the  winter  than  during  the  summer.  The  transition  between 
winter  to  summer  was  determined  by  the  change  in  the  sign  of  the  ECjC  and  occurred 
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rapidly  (<  1  month).  During  the  transitional  periods,  designated  here  as  the  spring  and 
fall,  the  eontribution  of  the  first  mode  to  the  mean  field  was  minimal. 

The  dominant  £[  wind  stress  pattern  was  centered  about  50°N,  165°W  and,  when 
combined  with  PCI ,  was  cyclonic  from  early-October  to  mid-March.  The  strongest 
annual  wind  stress  0.1  Pa)  were  southwestward  which  occurred  in  January  between 
30°-40°N  around  the  dateline.  These  were  comparable  to  magnitudes  of  the  mean  field 
found  in  this  region  (Figure  2)  however,  the  direction  of  the  stress  was  different.  The 
mean  wind  stress  field  showed  maximum  northwestward  wind  stress  (~  0.12  Pa)  located 
between  40°-50°N  east  of  the  dateline  that  extended  into  the  Gulf  of  Alaska.  The  mean 
field  wind  stress  in  this  region  was  part  of  a  subtropical  anticyclonic  circulation  whereas 
the  Ml  wind  stress  in  this  region  was  part  of  a  cyclonic  circulation  pattern  in  the  winter 
(October  -  March).  During  summer  (April  -  September)  the  dominant  wind  stress 
pattern  reversed  and  became  anticyclonic. 

Separate  and  less  dominant  features  were  also  visible  in  El  along  the  boundaries. 
In  the  Western  Pacific,  equatorward  wind  stress  associated  with  the  monsoons  occurred 
in  the  winter.  Specifically,  El  wind  stress  in  the  East  China  Sea  was  southward  while 
that  in  the  South  China  Sea  was  southwestward.  The  magnitude  (-0.08  Pa)  and  relative 
direction  during  the  winter  were  comparable  with  that  found  in  the  long-term  mean  wind 
stress  field.  In  the  Eastern  Pacific,  poleward  wind  stress  along  the  California-Oregon- 
Washington  coast  occurred  during  the  winter  and  reversed  equatorward  (upwelling 
favorable)  during  the  summer.  Also  in  the  Eastern  Pacific  there  was  a  region  of 
anticyclonic  circulation  centered  at  22.5°N,  120°W  west  of  the  Baja  California  Peninsula. 

The  Ekman  pumping  pattern  observed  in  Ey  (Eigure  9b)  appeared  related  to  the 
wind  stress  observed  in  El  (Eigure  9a).  In  fact,  straightforward  fluid  dynamical 
calculations  of  Ekman  pumping  from  El  wind  stress  gave  virtually  the  same 
geographical  pumping  pattern.  In  describing  E^  pumping,  it  was  observed  that  open- 
ocean  upwelling  (downwelling)  occurred  from  mid-October  to  March  over  much  of  the 
basin  north  (south)  of  the  zero  isotach  which  lies  roughly  along  30°N.  This  pattern 
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reversed  from  April  to  mid-October.  North  of  30°N,  the  interior  basin  of  the  mid¬ 
latitudes  and  the  Gulf  of  Alaska  regions  contoured  by  £'"=100  had  a  maximum  of  5.9 
cm  day  upwelling  in  January  and  5.1  cm  day  downwelling  in  August.  Pumping 
south  of  30°N  was  relatively  quiescent  in  the  mid-latitudes,  but  pumping  maxima  were 
found  in  the  tropics  along  15°N  and  west  of  the  dateline.  These  regions  showed 
maximum  downwelling  (~1 1.8  cm  day  in  January  and  maximum  upwelling  (-10.2  cm 
day  ”^)  in  August. 

The  pumping  pattern  along  the  western  boundary  consisted  of  intermixed  centers 
of  oppositely  phased  pumping.  For  instance,  along  the  western  boundary,  maximum 
(£*>200)  pumping  regions  occur  around  the  Kamchatka  Peninsula,  the  Japan  Sea,  and 
the  South  China  Sea.  Most  of  the  Asian  coast  north  of  30°N  was  an  upwelling  region 
during  the  winter  while  the  western  boundary  of  the  Bering  Sea,  east  of  the  Kamchatka 
Peninsula,  subtropical  and  tropical  coastal  regions  were  downwelling  during  this  time. 
The  pattern  was  reversed  in  summer. 

£"  pumping  along  the  eastern  boundary  had  a  less  complex  structure  than  that 
found  along  the  western  boundary.  The  zero  isotach  separated  nearshore  pumping  from 
oppositely-signed  interior  ocean  pumping  along  the  eastern  boundary  north  of  30°N.  A 
narrow  strip  of  nearshore  pumping  began  at  the  center  of  the  Gulf  of  Alaska  (140°W)  and 
continued  uninterrupted  southward  along  the  North  American  coast  to  the  southern  Baja 
Peninsula.  This  entire  boundary  pumping  region  was  in  phase  with  the  subtropical  open 
ocean  pumping;  downwelling  in  the  winter  and  upwelling  in  the  summer.  A  concentric 
region  of  intensified  pumping  was  found  within  this  region  and  was  located  from  the 
Columbia  River  southward  to  the  San  Francisco  Bay.  This  concentric  region  was  an  area 
of  intense  downwelling  during  the  winter  and  intense  upwelling  during  the  summer.  The 
maximum  was  located  at  about  42°N  (£"  =  300).  Maximum  downwelling  in  this  small 
region  occurred  in  January  at  a  rate  of  17.7  cm  day'^  and  maximum  upwelling  occurred  in 
August  at  a  rate  of  15.3  cm  day'\ 

The  most  striking  difference  in  comparison  with  the  long-term  mean  field  (Figure 
2b)  was  the  slope  of  the  zero  isotach  across  the  North  Pacific.  £"  had  a  zero  isotach  that 
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was  relatively  zonal  eentered  about  30°N.  The  long-term  mean  field,  however,  showed  a 
zero  isotaeh  that  began  in  the  west  about  30°N  and  sloped  northward  aeross  the  Paeifie  to 
the  Gulf  of  Alaska  (-'50°N)  before  turning  south  and  eneompassing  the  eastern  boundary. 
The  western  boundary  pumping  pattern  observed  in  the  long-term  mean  was  eonsistent 
with  the  winter-time  pattern  observed  in  with  regards  to  the  loeation,  approximate 
magnitude  and  phase  of  pumping  maxima.  For  example,  pumping  maxima  were 
observed  in  both  the  long-term  mean  and  in  the  Sea  of  Okhotsk  and  east  of  the 
Kamchatka  Peninsula,  in  the  Sea  of  Japan  and  eastward  of  Japan,  in  the  Yellow  Sea  and 
the  East  and  South  China  Seas.  Additionally,  the  position  of  the  zero  isotaeh  in  the 
Bering  sea  was  well  correlated  between  the  two  fields. 

Eastern  boundary  comparisons  between  the  long-term  mean  and  E^were  not  as 
well  correlated.  Most  of  the  midlatitude  and  subtropical  eastern  boundary  was  out  of 
phase  with  the  long-term  mean  field.  Additionally,  the  pumping  maxima  were  not  found 
in  the  same  position  nor  were  they  the  same  geographical  dimension.  The  concentric 
maxima  was  displaced  slightly  to  the  north  of  the  long-term  mean  maxima  found  along 
the  eastern  boundary.  Eong-term  mean  upwelling  along  the  eastern  boundary  extended 
from  northern  California  to  central  Mexico  and  was  more  intense  than  that  found  in  E^ 
pumping  at  any  time  of  the  year.  Still,  a  unique  nearshore-offshore  dipolar  upwelling- 
downwelling  feature  was  observed  in  both  fields  along  the  eastern  boundary  off  the 
California  coast,  but  in  slightly  different  locations.  Additionally,  pumping  activity 
around  the  Hawaiian  Islands  was  not  as  robust  in  E^  as  in  the  long-term  mean  pumping 
field. 

Ey  (Eigure  9c)  consisted  of  two  gyres  located  at  the  western  boundary.  The 
strongest  was  cyclonic  and  located  at  45°N.  This  gyre  was  centered  on  the  division 
between  the  subpolar  and  subtropical  gyre  as  defined  by  the  mean  field  (Eigure  2c).  The 
weaker  gyre  was  located  immediately  to  the  south  at  20°N  within  the  subtropical  gyre 
and  was  anticyclonic.  Maximum  amplitude  (0.05)  occurred  in  January  and  the  minimum 
in  September  (-0.04).  In  January,  would  increase  the  Kuroshio  transport  by  5.7  Sv 
and  in  September  decrease  it  by  3.9  Sv. 
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Figure  9.  Eigenvector  ( )  Plots,  {a)  El .  Scaling  Vector  is  0.25  Pa.  (b)  E".  Contour 
Interval  is  100  cm  Day’' .  {c)  .  Solid  Lines  Indicate  Negative  Values.  Arrowheads 

Indicate  Direction  of  Transport.  Zero  is  Denoted  by  the  Bold  Line. 

C.  EOF  MODE  2 

(Figures  10a, b  and  11a),  M'^  (Figures  10c  and  11b)  and  (Figures  lOd 
and  lie)  accounted  for  14.9%,  9.9%  and  24.5%  of  the  total  variance,  respectively. 
Whereas  the  first  principal  components  of  wind  stress,  Ekman  pumping  and  Sverdrup 
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transport  all  possessed  a  highly  eoherent  annual  eyele,  of  the  seeond  prineipal 
eomponents,  only  Ekman  pumping  {PC 2)  had  a  predominantly  annual  eyele.  Coherenee 
funetion  estimates  for  the  entire  monthly  averaged  time  series  (1948-1999)  for  PCI  and 
PC2  showed  the  two  time  series  to  be  highly  eoherent  at  the  semi-annual  eyele  (Figure 
8).  PCI ,  PC2  and  PC2  were  all  at  their  lowest  amplitude  in  November.  PCI  and 
PC2  had  larger  standard  deviations  during  winter  than  summer  while  the  standard 
deviation  of  PC 2  was  relatively  eonsistent  throughout  the  year. 


Figure  10.  Coherenee  Function  Fstimate  for  the  Entire  Monthly  Averaged  Time  Series 
(1948-1999)  of  Re(pC2)  and  PC^  (Solid  Fine),  Re(pCj)  and  PC^  (Broken  Fine), 

PC2  and  PC2  (Broken  Dotted  Fine). 

The  circulation  pattern  shown  for  El  (Figure  11a)  was  dominated  by  a  basin 
wide  circulation  centered  at  ~38°N/180°.  Cyclonic  circulation  occurred  in  winter 
(January  -  February)  and  summer  (May  -  September)  when  PCl  was  positive.  During 
this  time  poleward  circulation  occurred  in  the  Bering  Straits  and  the  South  and  East 
China  Seas.  Anticyclonic  circulation  occurred  during  the  spring  (March  through  April) 
and  the  fall  (September  to  mid-December)  while  there  was  equatorward  circulation  in  the 
Bering  Straits  and  the  South  and  East  China  Seas.  Wind  stress  was  greater  in  the 
subpolar  region  than  in  the  subtropical  region.  The  strongest  wind  stress  was  found 
between  45°-52°N.  Seasonally,  the  strongest  Ml  wind  stress  was  in  November  when 
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PCI  amplitudes  approach  -1.0.  The  negative  phase  of  Ml  was  stronger  in  fall  than 
spring  while  the  positive  phase  was  stronger  in  the  winter.  Although  these  seasonal 
patterns  are  clear  in  the  rotation  angle  and  magnitude  plots  (Figure  10a),  the  winter  and 
spring  transition  periods  are  within  one  standard  deviation  of  the  zero  line  so  that  the 
probability  that  the  second  mode  will  reverse  in  a  given  month  was  high.  Periods  of  time 
that  Ml  was  statistically  significant  (i.e.,  PCI  was  more  than  one  standard  deviation 
from  the  zero  line)  were  during  June  to  August  and  October  to  mid-November. 

The  most  striking  feature  of  (Figure  11b)  was  that  the  subpolar  and 
subtropical  pumping  regions  were  out  of  phase  with  midlatitude  and  eastern  boundary 
regions.  M^  midlatitude  pumping  covered  most  of  the  interior  North  Pacific  and  was 
distinctly  separated  from  the  eastern  boundary  pumping  region.  PC"  (Figure  10b)  was 
positive  for  the  first  seven  months  of  the  year  and  negative  for  the  last  five  months  of  the 
year.  The  maximum  amplitude  was  in  June  and  the  minimum  amplitude  was  in  October 
and  November.  The  resulting  M"  Ekman  pumping  was  strongest  during  the  fall  at 
which  time  it  was  twice  as  strong  as  during  the  summer.  Subpolar  pumping  maxima  were 
located  in  the  Sea  of  Okhotsk  and  the  Gulf  of  Alaska.  In  these  regions,  M"  pumping 
during  the  first  seven  months  of  the  year  was  downwelling  with  a  maximum  of  4.2  cm 
day"'  in  June  and  the  last  5  months  of  the  year  was  upwelling  with  maxima  during 
October  and  November  averaging  of  6.5  cm  day"'.  A  tropical  maxima  located  in  the 
western  Pacific  had  a  similar  phase,  but  twice  the  magnitude  (8.4  cm  day'^  in  June  and  13 
cm  day"'  in  October/November).  The  subpolar  and  tropical  regions  were  coincident  with 
active  pumping  regions  observed  in  P"  (Figure  9a)  and  the  mean  field  (Figure  2b). 

M"  pumping  along  the  eastern  boundary  and  in  regions  associated  with  the 
midlatitude  open  ocean  pumping  were  areas  of  upwelling  during  the  first  seven  months  of 
the  year  and  downwelling  during  the  last  five  months.  The  entire  eastern  boundary  from 
the  Queen  Charlotte  Islands  to  15°N  was  delineated  by  the  zero  isotach.  The  zero  isotach 
of  P"  (Figure  9b)  and  P"  was  in  the  same  position  along  the  eastern  boundary  from  the 
Queen  Charlotte  Islands  southward  to  the  San  Francisco  Bay  (Figure  1  lb).  Along  the 

eastern  boundary,  M"  and  M^  pumping  combined  to  enhance  mean  upwelling  by  as 
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much  as  15  cm  day'*  from  April  to  August  and  mean  downwelling  by  as  much  as  12  cm 
day'*  from  November  to  mid-December. 

E2  (Figure  19c)  consisted  of  two  gyres  located  at  the  western  boundary.  The 
strongest  was  anticyclonic  and  loeated  at  50°N.  Immediately  to  the  south,  a  eyclonie 
gyre  was  loeated  at  35°N.  The  northern  gyre  was  loeated  in  the  subpolar  gyre  and  the 
southern  gyre  in  the  subtropical  gyre  as  defined  by  the  mean  field  (Figure  2e).  PC^  had 
semiannual  variability  and  was  positive  from  mid-December  to  February  and  May  to 
mid-September.  The  standard  deviation  exeeeded  the  amplitude  in  all  months  exeept 
during  the  summer  months  of  June,  July  and  August.  During  this  time,  would 
inerease  the  Kuroshio  transport  by  2.2  Sv. 

D,  EOF  MODE  3 

The  third  mode  had  more  complex  spatial  variability  than  the  first  and  second 
mode.  Ml  (Figures  12a,  b  and  13a),  (Figures  12c  and  13b)  and  (Figures  12d 
and  13e)  account  for  8.5%,  5.1%  and  9.3  %  of  the  total  variance,  respectively.  The 
circulation  pattern  shown  for  El  (Figure  13a)  occurred  from  November  to  Mareh  as 
shown  by  the  rotation  angle  in  Figure  12a.  The  winter  circulation  pattern  was  strongest 
in  the  subpolar  region  where  flow  eonsisted  of  three  cireulation  centers,  antieyelonie  over 
the  Aleutian  Islands  and  cyelonic  over  Kamchatka  and  British  Columbia.  The  resulting 
eireulation  was  poleward  along  the  dateline  and  equatorward  along  the  coasts  of  Asia  and 
North  America.  In  the  subtropies,  the  strongest  eireulation  was  assoeiated  with  the 
monsoons  in  the  Western  Paeific  where  westward  flow  occurred  in  winter  but  weak 
eyelonie  flow  was  centered  over  Hawaii  and  weak  antieyelonie  flow  to  the  west  of  Baja 
California.  This  circulation  pattern  reverses  in  summer  (April  through  September)  so  that 
in  subpolar  regions,  equatorward  winds  oeeurred  along  the  dateline  and  poleward  winds 
along  the  eontinental  margins  and  in  the  Gulf  of  Alaska.  The  wind  stress  was  strongest 
during  the  period  November  through  January  and  in  May  and  July,  when  amplitudes 
(Figure  12a)  approaeh  0.4.  Although  the  seasonal  patterns  are  clear  in  the  rotation  angle 
(Figure  12a),  the  range  of  plus  or  minus  one  standard  deviations  erosses  the  zero  line  in 
all  months  (Figure  12b)  so  that  the  third  mode  will  reverse  in  a  given  month  was  high. 
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Figure  11.  Seasonal  Time  Series  of  .  (a)  rotationangleyPCl-)  is  the  Solid  Line  and 
magnitude{PClj^ )  is  the  Broken  Line  (b)  Re(F’C2^ )  is  the  Solid  Line  and  in  lm(F’C2^ ) 
the  Broken  Line,  (e)  (d)  PC'lL .  Shaded  Region  is  ±  Icr . 


Figure  12.  Eigenvector  ( )  Plots.  {&)  El .  Scaling  Vector  is  0.25  Pa.  (h)  .  Contour 

Interval  is  100  cm  Day"' .  (c)E'^.  Solid  Lines  Indicate  Negative  Values.  Arrowheads 
Indicate  Direction  of  Transport.  Zero  is  Denoted  by  the  Bold  Line. 

At  first  glance,  E^  (Ligure  13b)  appeared  related  to  El .  E^  (Ligure  13b)  was 
dominated  by  a  large  dipole  along  the  dateline  with  downwelbng  north  of  45  °N  centered 
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on  the  Aleutian  Islands  and  up  welling  to  the  south  centered  at  about  38°N.  Along  the 
eastern  and  western  boundaries,  this  pattern  reversed  with  downwelling  in  the  subpolar 
region  and  upwelling  in  the  subtropics  but  the  Western  Pacific  was  more  complex  and  the 
downwelling  area  along  28°N  extended  eastward,  almost  reaching  the  dateline.  A  small 
dipole  also  existed  in  the  region  of  the  Hawaiian  Islands;  however  the  variability  of  PC'^ 
(Figure  12b)  was  semiannual  (compared  to  the  annual  variability  of  the  rotation  angle  for 
PCI )  and  was  positive  in  January  and  February  and  from  July  through  October.  The 
maximum  amplitude  (0.04)  occurred  in  January  and  minimum  amplitude  in  April  (-0.04). 
As  for  PCI ,  the  standard  deviation  of  PC^  (Figure  12c)  was  equal  to  or  exceeded  the 
amplitude  in  each  month,  except  in  August. 

(Figure  13c)  consisted  of  four  gyres.  The  strongest  was  cyclonic  and  located 
at  the  western  boundary  at  56°N.  Immediately  to  the  south,  an  anticyclonic  gyre  was 
centered  at  48°N,  170°W.  Both  of  the  northern  gyres  were  located  within  the  subpolar 
gyre  as  defined  by  the  mean  field  (Figure  2c).  Two  weaker  gyres  occurred  farther  south 
within  the  subtropical  gyre,  cyclonic  at  34°N  and  anticyclonic  at  20°N.  PC^  (Figure 

12c)  had  annual  variability  and  was  positive  from  October  to  March  and  negative  the  rest 
of  the  time.  Maximum  amplitude  (0.025)  occurred  in  February  and  minimum  amplitude 
(-0.04)  in  July.  The  standard  deviation  of  PC'^  exceeded  the  amplitude  in  all  months 
except  July.  In  July,  M'^  would  increase  the  Kuroshio  transport  by  1 .2  Sv. 
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Figure  13.  Seasonal  Time  Series  of  .  (a)  rotationangle{PCl-)  is  the  Solid  Line  and 
magnitude{PClj^ )  is  the  Broken  Line  (b)  Re(F’C3^ )  is  the  Solid  Line  and  lm(F’C3'^^ )  is 
the  Broken  Line,  (c)  (d)  PC'^:^  .  Shaded  Region  is  ±  Icr . 
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Figure  14.  Eigenvector  ( )  Plots,  {a)  El .  Scaling  Vector  is  0.25  Pa.  (h)  E'^ .  Contour 
Interval  is  100  cm  Day’' .  (c)E'3^.  Solid  Lines  Indicate  Negative  Values.  Arrowheads 
Indicate  Direction  of  Transport.  Zero  is  Denoted  by  the  Bold  Line. 


IV.  DISCUSSION 


This  section  will  demonstrate  that  the  annual  variability  can  reasonably  be 
described  using  just  the  first  three  EOF  modes.  The  benefit  in  using  a  minimal  number  of 
modes  to  explain  annual  variability  is  that  the  complex  overlapping  patterns  modulated  in 
time  and  space  are  reduced  to  a  fewer  number  of  features  (from  the  eigenfield  spatial 
plots)  with  a  seasonal  amplitude  (from  the  principal  components). 

contains  the  largest  amount  of  total  variance  of  the  annual  cycle. 
had  an  annual  cycle  that  can  be  described  as  two  seasons  (winter  and  summer)  that  last 
for  five  months  each  separated  by  two  rapid  (approximately  one  month)  transition  periods 
(spring  and  fall).  North  Pacific  (NP)  winter  occurred  during  the  positive  phase  of 
(November  thru  March)  and  summer  during  the  negative  phase  (May  thru 
September).  The  transitional  periods  occurred  about  the  month  related  to  the  zero  line 
crossings,  April  and  October  for  the  spring  and  fall,  respectively.  Therefore,  seasons 
were  designated  based  on  amplitudes  and  zero  line  crossings. 

Individual  months  were  selected  to  represent  typical  seasonal  circulations  rather 
than  averaging  several  months  within  a  season.  Selecting  individual  months  provided  a 
truer  representation  of  actual  seasonal  circulation  patterns  and  magnitudes  at  smaller 
scales  which  may  have  been  smoothed  out  through  averaging.  The  selection  criteria  was 
based  on  the  minimum  and  maximum  amplitudes  of  and  where  the  zero  line 

crossings  occurred.  Maxima  PCl’'^’'^  amplitudes  occurred  in  January,  which  was 
selected  as  the  month  representative  of  winter.  Although  the  minima  for  PCI ,  PC^  ,  and 
PC^  occurred  at  different  times  of  the  year  (July,  August  and  September,  respectively), 
August  was  selected  as  the  month  representative  of  summer.  The  rationale  behind 
choosing  August,  other  than  it  was  the  median  month  among  the  three  circulations,  was 
that  the  August  value  was  within  one  standard  deviation  of  the  true  minima  for  each 
circulation,  and  therefore  only  magnitude,  not  phase,  would  be  slightly  affected  and 
remain  within  statistical  tolerance.  The  same  logic  was  used  for  choosing  April  and 
October  as  the  spring  and  fall  transition  months,  respectively.  Monthly  means 
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representative  of  the  winter  season  (a),  spring  transition  (b),  summer  season  (e)  and  fall 
transition  (d)  are  shown  for  wind  stress,  Ekman  pumping  and  Sverdrup  transport  in 
Figures  15,  19  and  23,  respeetively. 

This  diseussion  eontains  three  seetions,  a.  Wind  stress,  b.  Ekman  pumping  and  e. 
Sverdrup  transport.  Eaeh  seetion  will  eontain  a  diseussion  of  the  annual  eyele  using  a 
monthly  mean  representative  of  the  two  seasons  and  two  transition  periods,  quantify  how 
well  aeeount  for  the  seasonal  variability  and  examine  a  regional  example  to 

demonstrate  how  the  PCs  replieate  the  annual  eyele.  Easily,  the  interannual  variability  of 
the  PCs  is  eompared  to  environmental  indiees. 

A,  WIND  STRESS 

Monthly  means  for  the  winter  season  (a),  spring  transition  (b),  summer  season  (e) 
and  fall  transition  (d)  for  wind  stress  are  shown  in  Figure  15.  The  annual  eyele  of  North 
Paeifie  wind  stress  is  dominated  by  variations  in  intensity  and  meridional  movement  of 
the  subpolar  Aleutian  low  and  subtropieal  North  Paeifie  high.  The  Aleutian  low 
generates  eyelonie  wind  stress  and  dominates  North  Paeifie  eireulation  in  winter.  The 
North  Paeifie  subtropieal  high  generates  anticyelonie  wind  stress  that  dominates  the 
summer  eireulation.  Eaeh  eireulation  system  has  its  own  annual  eyele  in  the  North 
Paeifie  and  is  eonsidered  separately  below. 

The  strongest  monthly  mean  eyelonie  wind  stress  oeeurs  during  the  winter  in  the 
subpolar  region  and  is  eentered  about  50°N,  170°E  (Figure  15a).  Cyelonie  wind  stress 
oeeurs  over  the  North  Paeifie  north  of  27°N  from  the  eastern  to  the  western  boundaries. 
This  wind  stress  eireulation  pattern  diminishes  in  magnitude  and  geographieal  extent 
throughout  the  year,  beeoming  eentered  over  Alaska  during  the  spring  transition  (Figure 
15b),  and  disappearing  eompletely  in  the  summer  (Figure  15e).  It  and  re-develops  over 
Alaska,  east  of  its  winter  position,  during  the  fall  transition  (Figure  15d). 

The  annual  eyele  of  the  subtropieal  antieyelonie  wind  stress  shows  moderate  wind 
stress  eentered  about  30°N,  130°W  in  the  winter  (Figure  15a).  This  eireulation  beeomes 
eentered  at  about  32°N,  170°W  and  expands  basin- wide  during  the  spring  (Figure  15b), 
intensifying  equatorward  wind  stress  along  the  eastern  boundary  in  the  California  Current 
System  (CCS).  During  summer  the  eireulation  pattern  beeomes  eentered  about  40°N, 
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145°W  (Figure  15c).  Strong  equatorward  windstress  continues  in  the  CCS;  however, 
windstress  in  the  western  Pacific  subsides  and  becomes  poleward.  During  the  fall 
transition,  a  ridge  develops  connecting  eastern  and  western  anticyclonic  centers  about 
35°N,  140°W  and  35°N,  175°E  (Figure  15d)  which  create  westward  stress  in  the 
subtropics  and  eastward  stress  in  the  subpolar  regions. 

The  wind  stress  patterns  are  similar  to  each  other  during  both  spring  and  fall 
(Figures  15b  and  15d).  There  is  strong  anticyclonic  circulation  over  most  of  the  North 
Pacific  basin  centered  about  30°N  and  weak  cyclonic  circulation  centered  over  Alaska. 
Eastward  wind  stress  is  found  in  the  subpolar  regions  and  westward  wind  stress  is  found 
in  the  subtropical  regions.  A  nearly  zonal  band  of  strong  eastward  wind  stress  centered 
about  45°N  is  observed  when  cyclonic  and  anticyclonic  circulation  centers  exist 
concurrently  in  the  eastern  North  Pacific.  In  the  eastern  North  Pacific,  the  wind  stress 
gradient  is  strongest  between  35°-55°N  at  140°W  between  the  centers  of  these  two 
different  circulation  systems. 

At  smaller  scales  there  are  regional  differences  between  the  spring  and  fall.  Most 
notably,  anticyclonic  wind  stress  circulation  during  spring  (Eigure  15b)  appears  as  a 
single  coherent  pattern  that  is  centered  about  32°N,  170°W.  Although  the  fall  (Eigure 
15d)  also  exhibits  a  large  scale  anticyclonic  mean  circulation  across  most  of  the  basin,  it 
appears  to  be  a  composite  of  two  separate  anticyclonic  circulation  systems  centered  at 
35°N,  170°E  and  35°N,  140°W.  Additionally,  in  the  western  North  Pacific,  strong 
northeastward  wind  stress  is  centered  around  40°N  during  the  spring  and  during  the  fall 
eastward  windstress  is  centered  around  50°N  and  across  the  Sea  of  Okhotsk. 
Equatorward  wind  stress  along  the  eastern  boundary  is  more  intense  during  the  spring 
than  during  the  fall. 

Winter  and  summer  are  distinctly  different.  Winter  has  two  distinct  circulation 
centers  present  and  strong  southward  monsoonal  wind  stress  while  summer  shows  near 
uniform  anticyclonic  circulation  over  most  of  the  North  Pacific  and  northward  monsoonal 
wind  stress.  The  annual  cycle  of  the  windstress  in  the  monsoon  region  is  discussed 
below. 
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Figure  15.  Monthly  Mean  Wind  Stress  (Pa).  Sealing  Veetor  is  0.25  Pa. 
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Figure  16.  Ml^  3  +  X  for  the  Months  Indicated.  Scaling  Vector  is  0.25  Pa. 


Figure  17.  Ml 2  3  Residuals  for  the  Months  Indieated.  Sealing  Veetor  is  0.25  Pa. 
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Maximum  eastward  stress  occurs  in  the  western  and  central  Pacific  between  30°- 
45°N  during  the  winter  (Figure  15d),  averaging  0.18  Pa.  Maximum  equatorward  wind 
stress  in  the  CCS  occurs  during  spring,  averaging  0.15  Pa.  (Figure  15b).  The  weakest 
basin-wide  mean  scalar  wind  stress  field  magnitude  occurs  during  the  summer  (0.03  Pa) 
(Figure  15c)  and  the  strongest  during  the  winter  (0.07  Pa)  (Figure  15a). 

The  sums  of  Ml 2  3  +  X  (annual  mean)  for  January,  April,  August,  and  October 
are  shown  in  Figure  16.  The  difference  between  the  vector  components  of  the  monthly 
mean  and  Ml 22  +  X,  termed  Ml 22  residuals,  is  shown  in  Figure  17  and  was  used  as  a 

quantitative  measure  of  replicability  of  seasonal  wind  stress.  Upon  visual  inspection  of 
Figures  15  and  16,  the  three  modes  accounting  for  57.9%  of  the  total  variance  fairly 
replicates  the  main  physical  features  discussed  above.  Spring  (Figure  17b)  and  fall 
(Figure  17d)  have  larger  residuals  than  winter  (Figure  17a)  and  summer  (Figure  17c). 
This  is  in  part  due  to  the  contribution  of  mode  1,  which  accounts  for  34.5%  of  the  total 
variance,  being  near  zero  at  these  times.  This  implies  the  higher  modes  may  account  for 
much  of  the  transitional  period  variability.  Regions  of  large  deviation  (in  parenthesis) 
occur  in  the  Bering  Sea  (0.04  Pa),  the  western  Pacific  (0.06  Pa),  the  Gulf  of  Alaska  (0.04 
Pa),  and  the  trade  wind  regions  (0.05  Pa)  during  the  spring  and  fall.  This  deviation  is 
accounted  for  in  the  higher  modes,  especially  Ml  and  Ml  (Appendix  D). 

Two  wind  stress  time  series  for  the  northwest  Pacific  are  shown  in  Figure  18. 
Geographic  boundaries  are  0°-50°N  and  100“-160°E  for  both  time  series.  Data  limits  are 
15°-50°N  and  100°-160°E  to  correspond  with  the  established  EOE  analysis  limits.  Wind 
stress  is  computed  on  a  2.5'’  grid  and  plotted  on  a  5°  grid  for  resolution  of  the  wind  stress 
vectors. 

The  Aleutian  low  is  fully  developed  during  the  winter,  generating  large  cyclonic 
wind  stress  over  much  the  western  North  Pacific.  The  Asian  winter  high  extends  over 
much  of  the  ocean  and  produces  a  wind  reversal  over  the  East  and  South  China  Seas  and 
the  region  east  of  the  Philippines;  these  regions  thus  experience  monsoonal  climate,  with 
Northeast  Monsoon  during  winter  (December  -  March)  and  Southwest  Monsoon  during 
summer  (June  -  September). 
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PCI  +  X  replicates  the  monthly  means  throughout  the  annual  cycle  well  even 
though  PCI  accounts  for  only  34.5%  of  the  total  variance  of  the  wind  stress.  During  the 
periods  of  May-August,  PCI  +  X  does  not  have  the  same  direction  as  the  monthly  mean 
although  the  magnitude  of  the  wind  stress  is  similar.  This  is  also  reflected  in  the  residual 
plots  (Figure  17)  and  may  indicate  that  higher  modes  are  important  during  this  time.  Fall 
and  winter  PCI  +X  wind  stress  show  similar  direction,  but  with  a  lesser  magnitude. 
This  is  also  apparent  in  the  residual  plots. 
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B,  EKMAN  PUMPING 

Since  wind  stress  on  the  oeean  surfaee  varies  in  time  and  spaee,  so  does  the 
Ekman  transport.  Diminishing  eyelonie  wind  stress  and  expanding  antieyelonie  wind 
stress  from  winter  to  summer  results  in  a  general  seasonal  northward  movement  of  the 
zero  isotaeh  (Figures  19a,  b,  e).  The  zero  isotaeh  is  at  its  southernmost  limit  (~32°N)  in 
the  western  and  eentral  Paeifie  during  winter  (Figure  19a)  and  its  northernmost  limit 
(~50°N)  during  summer  (Figure  19e).  The  mean  upwelling  value  for  the  region  north  of 
the  zero  isotaeh  is  greatest  (-'5.9  em  day'^)  in  the  winter  and  weakest  in  the  summer  (-1.4 
em  day'^). 

Sub-basin  seale  upwelling  maxima  in  the  Gulf  of  Alaska  and  near  the  Kamehatka 
peninsula  persist  throughout  the  year,  but  are  strongest  (>  15  em  day'^)  when  eyelonie 
wind  stress  generated  from  the  Aleutian  low  is  present  during  winter  (Figure  19a)  and  fall 
(Figure  19d).  A  distinetive  tilt  of  the  zero  isotaeh  oeeurs  in  the  Gulf  of  Alaska  during  the 
winter  when  the  Aleutian  Fow  is  eentered  west  of  the  dateline.  Weak  (<  5  em  day'^) 
small-seale  coastal  downwelling  also  occurs  off  Oregon,  Washington  and  British 
Columbia  during  these  seasons. 

A  small-seale  upwelling  feature  near  the  Hawaiian  Islands  (20°N,  160°W),  that  is 
also  apparent  in  other  elimatologies  (Trenberth  et  al.,  1990;  Chelton  et  al.,  1990; 
Hellerman  and  Rosenstein,  1983),  persists  throughout  the  year  and  is  due  to  ehanges  in 
both  wind  speed  and  direetion  (Figures  15a-d).  The  formation  of  a  dipolar  upwelling- 
downwelling  strueture  in  the  summer  (Figure  19e)  and  fall  (Figure  19d)  has  been  noted 
and  diseussed  in  other  studies  (Trenberth  et  al,  1990;  Harrison,  1989). 
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The  CCS  is  dominated  by  intense  equatorward  alongshore  (upwelling  favorable) 
wind  stress  during  the  spring  into  the  summer  and  fall  (Figures  15b-d).  Coastal 
upwelling  in  the  CCS  is  most  intense  (>  15  cm  day'^)  off  the  California  and  Baja  coasts 
during  the  spring  (Figure  19b).  The  active  upwelling  region  covers  the  largest  zonal 
range  in  summer  (Figure  19c),  extending  from  British  Columbia  to  the  Baja  Peninsula.  In 
summer,  areas  with  upwelling  rates  greater  than  15  cm  day'^  are  located  from  Oregon  to 
the  central  California. 

This  upwelling  favorable  wind  stress  is  associated  with  anticyclonic  curl  offshore 
of  the  coastal  upwelling  zone  in  the  CCS  (Figure  15).  The  pairing  of  coastal  cyclonic 
curl  and  offshore  anticyclonic  curl  results  in  the  development  of  a  dipolar  upwelling- 
downwelling  structure  (30°-45°N,  120°-135°W)  (Figure  20).  This  dipole  oscillates  north 
and  south  seasonally  about  35°N.  The  dipole  feature  intensifies  as  it  moves  northward 
through  the  spring  transition  and  summer.  The  dipole  is  strongest  and  at  its  northernmost 
limit  (~  45°N)  in  July  and  at  its  southernmost  limit  (~  32°N)  in  December  thru  February. 

The  offshore  downwelling  feature  tracks  the  northward  movement  of  the 
upwelling  maximum  and  intensifies  at  the  same  time.  The  upwelling  portion  of  the 
dipole  is  most  intense  in  the  months  of  June  and  July  (summer),  averaging  about  35  cm 
day'\  and  least  intense  from  December  through  February  (winter),  averaging  about  15  cm 
day''.  The  downwelling  portion  of  the  dipole  is  most  intense  during  June  and  July  also, 
averaging  25  cm  day''.  This  feature  averages  about  15  cm  day''  from  August  through 
March  and  20  cm  day''  during  April  and  May. 

The  implication  of  all  this  is  that  the  seasonal  evolution  of  coastal  upwelling 
relates  to  the  large-scale  evolution  of  the  wind  field  over  the  North  Pacific.  Also  the 
onset  of  spring  transition  is  due  to  the  variability  in  seasonal  transition  from  an  Aleutian 
low  to  a  North  Pacific  high  dominance. 
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Figure  20.  Monthly  Mean  Ekman  Pumping  or  Vertieal  Velocity  at  the  Base  of  the  Ekman 
Eayer  as  a  Result  of  Divergence  of  the  Ekman  Transport  Computed  for  the  Period  1948- 
1999.  The  Contour  Interval  is  5  cm  day'\  Zero  Isotach  is  Bold.  Upwelling  Regions  are 

Shaded.  Geographic  Eimits  As  Shown. 

The  sum  of  ^  (annual  mean)  for  January,  April,  August,  and  October 

are  shown  in  Figure  21.  The  difference  between  the  monthly  mean  and  ^ 

termed  Mj"  3  residuals,  is  shown  in  Figure  22  and  was  again  used  as  a  quantitative 

measure  of  replicability.  The  three  modes,  accounting  for  only  38.5%  of  the  total 
variance,  replicate  the  main  physical  features  discussed  above  to  within  5  cm  day'\  The 
CCS  dipole,  the  Hawaii  Island  dipole  and  active  upwelling  centers  are  well  replicated. 
The  winter  tilt  of  the  zero  isotach,  sloping  northeastward  east  of  155°W,  is  well  replicated 
also. 


43 


C.  SVERDRUP  TRANSPORT 

The  total  ocean  transport  based  on  the  Sverdrup  relation  (Sverdrup,  1947; 
Pedlosky,  1987)  is  used  as  a  measure  of  the  oceanic  response  to  atmospheric  forcing, 
rather  than  a  true  reflection  of  the  actual  currents.  Calculations  of  the  Sverdrup  transport 
are  meaningful  when  the  forcing  is  slowly  varying,  relative  to  the  internal  adjustment 
time  of  the  ocean,  i.e.  the  assumption  of  a  steady  state  is  valid.  The  reliability  of  this 
assumption  decreases  with  increasing  latitude.  Even  when  the  wind  stress  varies  with 
time,  it  must  be  assumed  that  the  ocean  is  always  in  near-equilibrium,  namely  the 
Sverdrup  balance  (Kutsuwada,  1988). 

Monthly  means  representative  of  the  winter  season  (a),  spring  transition  (b), 
summer  season  (c)  and  fall  transition  (d)  for  Sverdrup  transport  are  shown  in  Figure  23. 
The  calculated  transports  are  in  general  agreement  with  Trenberth  et  al.  (1990)  and 
Hellerman  and  Rosenstein  (1983)  (Table  3).  Differences  arise  due  to  the  different  drag 
coefficients  and  wind  data  sets  used. 

Seasonal  variations  in  the  wind  stress  change  the  strength  of  the  subpolar  and 
subtropical  gyres  substantially  (Trenberth  et  ah,  1990).  The  annual  mean  eastward 
transport  derived  from  NCEP  winds  is  80  Sverdrups  (Figure  2c).  The  gyres  are  strongest 
in  northern  hemisphere  winter  (Figure  23  a)  and  weakest  in  summer  (Figure  23  c). 
Seasonal  variations  range  from  128  Sv  in  January  to  around  56  Sv  in  August.  The 
circulation  of  the  North  Pacific  subpolar  gyre  decreases  from  57  to  33  Sv  (-42%)  and  the 
subtropical  gyre  decreases  from  71  to  23  Sv  (-68%)  accounting  for  most  of  the  decrease 
in  the  eastward  transport. 
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Winter 

NCEP  EC 

HR 

Summer 

NCEP  EC 

HR 

Annual 

NCEP  EC 

HR 

Northward  Transport 

Subpolar  gyre 

-57 

-75 

-55 

-33 

-25 

-10 

-40 

-50 

-35 

Subtropical  gyre 

71 

100 

100 

23 

40 

45 

40 

60 

60 

Eastward  Transport 

North  Pacific 

128 

175 

155 

56 

65 

55 

80 

no 

95 

Table  3.  North  Pacific  Transports  (Sv)  Values  based  on  NCEP  Wind  Stress  Forcing  from 
1948-1999  for  Winter  (February)  and  Summer  (September).  For  Comparison, 
Transports  from  Two  Other  Studies  for  Winter  (January)  and  Summer  (July)  are 
Shown.  Sources  are  Trenberth  Ft  Al.  (1990),  Derived  from  FCMWF  (FC)  and 
Hellerman  and  Rosenstein  (HR)  Wind  Climatologies. 

The  sum  of  ^  (annual  mean)  for  January,  April,  August,  and  October 

are  shown  in  Figure  24.  The  Mj"  3  residuals  are  not  shown  as  the  residual  values  are  less 

than  10  Sv  for  each  of  the  four  months.  In  other  words,  the  first  three  modes,  accounting 
for  only  74.6%  of  the  total  variance,  replicate  the  main  physical  features  discussed  above 
to  within  10  Sv. 
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While  the  Sverdrup  relation  does  not  replicate  the  western  boundary  current  (i.e. 
the  Kuroshio  current),  it  is  assumed  that  the  transport  of  the  mid-latitude  anticyclonic 
gyre  must  be  balanced  by  a  western  boundary  current.  The  circulation  inferred  from  the 
streamfunctions  in  the  subtropics  gives  rise  to  the  Kuroshio  as  the  return  flow.  The 
Kuroshio  transport  is  defined  as  the  total  eastward  flow  north  of  30°N  (Qiu  and  Joyce, 
1992).  This  transport  is  composed  of  a  barotropic  (Sverdrup)  and  baroclinic 
(geostrophic)  components.  The  Kuroshio  transport  has  a  large  seasonal  signal  (Sekine 
and  Kutsuwada,  1994),  with  maximum  baroclinic  transport  in  the  summer;  the  maximum 
Sverdrup  transport  calculated  from  the  winds  occurs  in  winter. 

Kutsuwada  (1988)  suggests  that  another  way  to  derive  the  Kuroshio  transport  is  to 
estimate  it  from  sea  level  data  at  two  stations  on  both  sides  of  a  strait  through  which  the 
Kuroshio  always  flows.  The  sea  level  variations  are  due  to  the  variations  both  of  the 
barotropic  transport  and  the  baroclinic  transport.  The  total  volume  transport  of  the 
Kuroshio  is  composed  of  transports  through  the  Osumi  and  Tokara  Straits  (Kawabe, 
1995). 

Figure  25  looks  at  the  time  variations  of  the  depth  averaged  transport  derived 
from  the  monthly  averaged  sea  level,  the  Sverdrup  transport  and  the  PC^  +  X  Sverdrup 
transport.  Sea  level  data  from  the  Joint  Archive  for  Sea  Level  (JASL)  at  the  University 
of  Hawaii  was  obtained  for  two  stations  on  opposite  sides  of  the  Tokara  Strait.  Monthly 
mean  values  from  1964  thru  1968  were  used  for  Naze,  Japan  (28°22.7’N,  129°29.9’E) 
and  Nishinoomote,  Japan  (30°43.9’N,  130°59.7’E).  A  mean/between  the  two  stations 
was  used  to  calculate  the  transport  at  a  depth  of  1000  meters.  The  Sverdrup  transport  and 
the  PC^  +  X  Sverdrup  transport  is  the  transport  between  the  two  grid  points  nearest  the 

tide  stations,  27.5°N,  130°E  and  30°N,  130°E.  The  Sverdrup  transport  is  derived  from 
NCEP  wind  observations  (1948-1999)  as  described  above 

There  is  some  discrepancy  between  the  barotropic  and  baroclinic  annual  cycle 
which  may  be  due  to  unmeasured  barotropic  changes  as  was  suggested  by  Hautala  et  al. 
(1994)  for  the  difference  between  the  annual  mean  Sverdrup  and  Kuroshio  transport  in 
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the  western  Paeifie.  The  Sverdrup  relation,  however,  aeeounts  for  a  good  portion  of  the 
seasonal  varianee  of  the  Kuroshio  transport,  giving  eredibiiity  to  the  results  shown  here. 


Figure  25.  Comparison  of  Depth  Averaged  Transport  in  the  Tokara  Strait  and  Sverdrup 
Transport  in  the  Region.  Seasonal  Variability  of  Depth  Averaged  Transport  (Thick  Black 
Line).  The  Sverdrup  Transport  (Thick  Grey  Line)  and  the  PC^  +  X  Sverdrup  Transport 
(Thin  Grey  Line)  is  the  Transport  between  the  Two  Grid  Points  Nearest  the  Tide  Stations, 
27.5°N,130°E  And  30°N,130°E.  The  Sverdrup  Transport  is  Derived  from  NCEP  Wind 
Observations  (1948-1999)  as  Described  Above 


D,  INTERANNUAL  VARIABILITY 

This  thesis  has  characterized  the  annual  variability  of  the  wind  stress,  Ekman 
pumping  and  Sverdrup  transport  in  the  North  Pacific.  The  EOE  analysis  may 
characterize  interannual  variability  of  these  circulation  fields  as  well.  The  PCs  were 
compared  with  an  atmospheric  climate  index  and  the  first  modes,  although  strongly 
correlated  with  each  other,  were  poorly  correlated  with  the  climate  index.  Eor 
illustration,  Eigure  26  shows  smoothed  time  series  of  PC^,  PC'^  and  PC 2  ■  Down  and  up 
arrows  indicate  the  time  of  the  strong  historic  El  Nino  and  Ea  Nina  events,  respectively 
(Wolter  and  Timlin,  1998). 


The  PCs  are  well  correlated  to  each  other  on  the  interannual  time  scale; 
'(^q^pc3>)=0•^6,  =  Strong  El  Nino  (Ea  Nina)  events  are 


consistently  characterized  by  increases  (decreases)  in  the  smoothed  time  series  of  all 
three  PCs  shown  in  Eigure  26.  The  implication  of  this  is  that  El  Nino  events  result  in  a 
stronger  cyclonic  wind  stress  around  the  Aleutian  Eow  (Eigure  11a),  and  a  weaker  North 
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Pacific  Current  (Figure  lie).  The  assoeiated  Ekman  pumping  will  be  greater  in  the 
eoastal  Gulf  of  Alaska  and  south  of  45°N,  and  weaker  in  the  California  Current  and 
Bering  Sea  (Figure  13b).  The  eneouraging  eorrelation  of  these  PCs  with  the  El  Nino- 
Southem  Oseillation  (ENSO)  signals  suggests  EOE  analysis  of  seasonal  values  of  these 
fields  will  yield  an  improved  understanding  of  the  response  of  the  North  Paeifie  Oeean  to 
ENSO  foreing. 


Eigure  26.  Smoothed  Time  Series  (36  point)  of  PC^  as  Indieated  Versus  Smoothed  Time 

Series  (36  point)  Monthly  N01(/-15)  Values  (Sehwing  et  al.,  2001)  (heavy  gray  line). 
NOl  is  an  Index  of  North  Paeifie  Climate  Variability. 


There  is  also  some  evidenee  of  eoherent  longer-term  variability  in  these  results. 
The  three  PCs  are  plotted  against  monthly  Northern  Oseillation  Index  (NOl)  values,  an 
environmental  index  of  interannual  to  deeadal  elimate  variability  in  the  North  Paeifie 
(Sehwing  et  ah,  2001).  The  PCs  follow  the  general  trend  of  the  NOl  and  explain  about 
one  quarter  to  one  third  of  the  interannual  variability  of  the  NOl  using  this  method  of 


EOE  analysis;  n 


{PCI,N0I) 


=  0.38,  n 


{pCi  ,NOl) 


=  0.25,  K 


(PCI'.VO/) 


=  0.36.  There  appears  to  be 


some  eoineidenee  between  the  PCs  and  the  NOl  on  deeadal  seales.  Eor  example,  all  series 
inerease  at  about  the  time  of  a  reported  regime  shift  in  1976  (Trenberth,  and  Hurrell, 
1994).  Euture  work  is  planned  using  a  similar  EOE  methodology  to  explore  interannual 
variability  by  season. 
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APPENDIX  A.  VARIANCE 


Percentage  of  total  variance  and  cumulative  variance  explained  by  the  first  10 

2  2  2 

EOF  modes  for  wind  stress  ( r  =  Pa  ),  Ekman  pumping  ( w  =cm  day'  )  and  Sverdrup 
transport  (i//  =  Sv^). 


Mode  r  w  ¥ 


Variance  Variance  Variance 


Percent 

Cumulative 

Percent 

Cumulative 

Percent 

Cumulative 

1 

34.5 

34.5 

23.5 

23.5 

45.6 

45.6 

2 

14.9 

49.4 

9.9 

33.4 

19.7 

65.3 

3 

8.5 

57.9 

5.1 

38.5 

9.3 

74.6 

4 

5.9 

63.8 

4.0 

42.5 

5.6 

80.2 

5 

5.1 

68.9 

2.9 

45.4 

3.2 

83.4 

6 

3.4 

72.3 

2.5 

47.9 

2.9 

86.3 

7 

2.4 

74.7 

2.2 

50.1 

1.7 

88.0 

8 

1.9 

76.6 

1.8 

51.9 

1.7 

89.7 

9 

1.3 

77.9 

1.8 

53.7 

1.4 

91.1 

10 

1.2 

79.1 

1.4 

55.1 

1.2 

92.3 
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APPENDIX  B.  LONG-TERM  MONTHLY  MEAN  PRINCIPAL 

COMPONENTS 


Long-term  (1948-1999)  monthly  mean  (m)  principal  components  (PC^-)  for 
wind  stress  (X  =  t),  Ekman  pumping  {X  =  w)  and  Sverdrup  Transport  {X  =  ^)  for  the 
first  (  PCy  ),  second  (  PCf  )  and  third  ( PC^  )  principal  components. 


Month 

PC^ 

PC^ 

PCf 

r 

w 

W 

T 

w 

W 

T 

w 

iP 

J 

1.375 

0.059 

0.52 

0.576 

0.012 

0.17 

0.324 

0.037 

0.25 

F 

1.178 

0.048 

0.48 

0.459 

0.025 

0.09 

0.172 

0.018 

0.25 

M 

0.508 

0.026 

0.25 

-0.017 

0.029 

-0.13 

0.148 

-0.021 

0.05 

A 

-0.256 

-0.001 

-0.04 

-0.306 

0.023 

-0.17 

-0.038 

-0.035 

-0.06 

M 

-0.542 

-0.018 

-0.07 

0.052 

0.031 

0.01 

-0.241 

-0.024 

-0.20 

J 

-0.811 

-0.033 

-0.09 

0.477 

0.042 

0.22 

-0.281 

-0.002 

-0.30 

J 

-1.135 

-0.047 

-0.27 

0.343 

0.020 

0.22 

-0.118 

0.006 

-0.38 

A 

-1.111 

-0.051 

-0.35 

0.279 

-0.000 

0.26 

-0.113 

0.013 

-0.13 

S 

-0.682 

-0.042 

-0.35 

0.071 

-0.028 

0.15 

-0.291 

0.019 

-0.04 

0 

-0.140 

-0.014 

-0.30 

-0.618 

-0.064 

-0.19 

-0.142 

0.004 

0.14 

N 

0.484 

0.021 

-0.07 

-0.998 

-0.066 

-0.41 

0.262 

-0.017 

0.22 

D 

1.130 

0.051 

0.30 

-0.318 

-0.026 

-0.22 

0.317 

0.002 

0.21 
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APPENDIX  C.  MAGNITUDE  AND  ROTATION  ANGLE 


Long-term  (1948-1999)  monthly  mean  magnitude  {PCl-)  (Amplitude)  and 
rotation  angle  {PCI-^  (Rotation)  for  the  first  {PCD,  second  {PCD  and  third  {PCD 
principal  components  of  wind  stress. 


Month 

PC 

r 

1 

PC 

r 

2 

PC 

•T 

Amnlitude 

Rotation 

Amplitude 

Rotation 

Amplitude 

Rotation 

J 

1.6537 

33.7229 

0.5815 

7.6212 

0.3344 

14.6895 

F 

1.3977 

32.5308 

0.4671 

10.9634 

0.1725 

355.6761 

M 

0.5946 

31.4013 

0.0302 

123.5166 

0.1684 

331.6285 

A 

0.3191 

216.6722 

0.3145 

193.6069 

0.0402 

160.4186 

M 

0.6227 

209.5626 

0.0581 

333.0765 

0.2830 

211.7337 

J 

0.9158 

207.6332 

0.4842 

10.0637 

0.3475 

216.0875 

J 

1.2784 

207.4295 

0.3443 

4.8415 

0.2019 

234.2814 

A 

1.2632 

208.4240 

0.2800 

4.4830 

0.1639 

226.6040 

S 

0.8140 

213.1499 

0.0923 

40.1402 

0.2972 

191.4542 

0 

0.2795 

239.9615 

0.6182 

182.1053 

0.2839 

119.8881 

N 

0.5008 

14.7813 

1.0138 

190.0247 

0.3804 

46.5604 

D 

1.3257 

31.4950 

0.3294 

195.2547 

0.3590 

27.8798 
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APPENDIX  D. 


Seasonal  time  series  of  the  PC  pairings,  referred  to  as  modes  and 

designated  .  The  prineipal  eomponents  (PC)  for  a  partieular  mode  (M)  and  the 

eorresponding  eigenveetors  (E)  are  eaeh  on  separate  pages.  Captions  at  the  bottom  of 
eaeh  page  referenee  the  mode. 

There  are  four  prineipal  eomponent  plots  per  mode:  (a)  rotation  angle{PC^^)  is 
the  solid  line  and  magnitude{PCl-)  is  the  broken  line  (b)  Re(E’C^)  is  the  solid  line  and 
Im(E’C^-)  is  the  broken  line,  (c)  PC^-{d)  PC^- .  Shaded  region  is  ±lcr. 

There  are  three  eigenveetor  plots  per  mode:  (a)  .  Sealing  veetor  is  0.25  Pa. 

(b)  E^ .  Contour  interval  is  100  em  day'\  (e)  E^ .  Solid  lines  indieate  negative  values. 
Arrowheads  indieate  direetion  of  transport. 
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Rotation  angle  (-) 


Seasonal  Time  Series  of  PC^- .  (a)  rotation  angleyPClj^j  is  the  Solid  Line  and 
magnitude{PCl- )  is  the  Broken  Line  (b)  Re(T’C4- )  is  the  Solid  Line  and  Im(T’C4- ) 
the  Broken  Line,  (e)  T’Cf-  (d)  PCI-.  Shaded  Region  is  ±lcr. 


Magnitude  (- 


Eigenvector  (£’f)  Plots,  (a)  .  Scaling  Vector  is  0.25  Pa.  (b)  EJ.  Contour  Interval  is 

lOOcmDay'^  {c)  .  Solid  Lines  Indicate  Negative  Values.  Arrowheads  Indicate 

Direction  of  Transport.  Zero  is  Bold. 
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IS 

le  Solid 
d  Reeic 


Eigenvector  (£’5'^)  Plots,  (a)  .  Scaling  Vector  is  0.25  Pa.  (h)  .  Contour  Interval  is 

100cmDay'\  (c)  .  Solid  Lines  Indicate  Negative  Values.  Arrowheads  Indicate 

Direction  of  Transport.  Zero  is  Bold. 
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Seasonal  Time  Series  of  PC^j^ .  (a)  rotation  angle{PClj^^  is  the  Solid  Line  and 
magnitude{PCljf^ )  is  the  Broken  Line  (b)  Re(T’Cg^ )  is  the  Solid  Line  and  \m{pCl^ ) 
the  Broken  Line,  (e)  T’CIL  (d)  PCI-.  Shaded  Region  is  ±lcr. 


100  140  ISO  140  100 


Eigenvector  (£’^)  Plots,  (a)  Eg  .  Scaling  Vector  is  0.25  Pa.  (b)E^.  Contour  Interval  is 
100cmDay'\  (c)  .  Solid  Lines  Indicate  Negative  Values.  Arrowheads  Indicate 

Direction  of  Transport.  Zero  is  Bold. 
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Rotation  angle  <— > 


Magnitude  (— 


Eigenvector  (£’^)  Plots,  (a)  .  Scaling  Vector  is  0.25  Pa.  (b)  E".  Contour  Interval  is 

100cmDay'\  (c)  Ef  .  Solid  Lines  Indicate  Negative  Values.  Arrowheads  Indicate 

Direction  of  Transport.  Zero  is  Bold. 


67 


Eigenvector  (£’g^)  Plots.  {^)  El .  Scaling  Vector  is  0.25  Pa.  (b)E’g‘^.  Contour  Interval  is 
100cmDay'\  {c)  E'^  .  Solid  Lines  Indicate  Negative  Values.  Arrowheads  Indicate 

Direction  of  Transport.  Zero  is  Bold. 


69 


Rotation  angle  <— > 


% 


Magnitude  (— 


100 


140 


ISO 


140 


100 


Eigenvector  Plots,  (a)  Eg  .  Scaling  Vector  is  0.25  Pa.  (h)  Eg  .  Contour  Interval  is 
100cmDay'\  (c)  Eg  .  Solid  Lines  Indicate  Negative  Values.  Arrowheads  Indicate 

Direction  of  Transport.  Zero  is  Bold. 
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Seasonal  Time  Series  of  PCf^-  .  (a)  rotation  angle  (t’Cjq-)  is  the  Solid  Line  and 
magnitude  [pCIq-  )  is  the  Broken  Line  (b)  Re(T’Cjo- )  is  the  Solid  Line  and  Im(T’Cjo- ) 
the  Broken  Line,  (e)  PCq- (d)  PCiq- .  Shaded  Region  is  ±lcr. 


100 


140 


ISO 


140 


100 


Eigenvector  ( 
is  100  cm  day''. 


)  Plots,  (a)  Ejo.  Scaling  Vector  is  0.25  Pa.  (b)  E".  Contour  Interval 
(c)  Ejq  .  Solid  Lines  Indicate  Negative  Values.  Arrowheads  Indicate 
Direction  of  Transport.  Zero  is  bold. 
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APPENDIX  E.  WIND  STRESS  MONTHLY  MEANS 


The  data  shown  in  this  appendix  are  monthly  means  of  wind  stress  (Pa  =  kg  m"' 
s' ).  The  values  were  derived  from  NCEP  daily  averaged  surface  winds  from  1948-1999. 
Monthly  mean  fields  are  computed  as  the  arithmetic  mean  of  each  month  over  the  entire 
time  series.  Scaling  vector  is  0.5  Pa. 
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APPENDIX  F.  EKMAN  PUMPING  MONTHLY  MEAN 


The  data  shown  in  this  appendix  are  monthly  means  of  Ekman  pumping  values. 
The  values  were  derived  from  monthly  averaged  wind  stress  values  from  the  NCEP 
subset  values  of  daily  wind  stress.  Climatology  fields  are  computed  as  the  arithmetic 
mean  of  each  month  over  the  entire  time  series.  The  zero  isotach  is  bold.  The  contour 
interval  is  5  cm"'  day.  Shaded  regions  indicate  positive  values  (upwelling). 
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APPENDIX  G.  SVERDRUP  TRANSPORT  MONTHLY  MEANS 


The  data  shown  in  this  appendix  are  monthly  means  of  Sverdrup  transport  values. 
The  values  were  derived  from  monthly  averaged  wind  stress  values  from  the  NCEP 
subset  values  of  daily  wind  stress.  Monthly  mean  fields  are  computed  as  the  arithmetic 
mean  of  each  month  over  the  entire  time  series.  The  zero  streamline  is  in  bold. 
Streamline  interval  is  10  Sv  (1  Sv  =  10^  m^  s"'). 
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APPENDIX  H.  COMPLEX  EMPIRICAL  ORTHOGONAL 
FUNCTION  (EOF)  MATLAB  PROGRAM 


87 


%  wind  Stress  Vector  Complex  EOF  calculation 
%10/27/99 

q, 

o 

%DEMEANING  AND  DETRENDING  U&V  INDIVIDUALLY  THEN  LEGLER 
callstart  %clear  all  figures,  variables 


q, 

o 


%load  u,v  component  monthly  matrices  (29,73,624) 

%52  years  of  monthly  averaged  u  (tau  x.mat)  and  v  (tau  y.mat) (624) 
%2 . 5  degree  grid  (29,73) 

load  tau_x.mat 
load  tau_y.mat 


q, 

o 

%create  empty  matrix 

new_tau_x  =  zeros ( [23  73  624]); 
new_tau_y  =  zeros ( [23  73  624]); 


%trim  data  from  0-70N  to  15-70N 
for  n=l : 624 

new_tau_x ( : , : , n)  =  tau_x(l:23,:,n); 
new_tau_y ( : ,  : , n)  =  tau_y ( 1 : 23 ,  : ,  n) ; 
n 

end 

q, 

o 

%DETREND,  DEMEAN  U,V  COMPONENTS 
%create  empty  matrix 

transpose_tau_x  =  zeros ([73  23  624]); 
transpose_tau_y  =  zeros ([73  23  624]); 
reshape_tau_x  =  zeros ([624  1679]); 

reshape_tau_y  =  zeros ([624  1679]); 


for  n  =  1:624 

transpose_tau_x ( : , : 

%TRANSPOSE 
transpose_tau_y ( : , : 

%TRANSPOSE 
reshape_tau_x (n, : ) 

% RE SHAPE 

reshape_tau_y (n, : ) 
%RESHAPE 
end 

q, 

o 

%create  empty  matrix 
dmn  reshape  x 
dmn  reshape  y 
dt  dmn  rs  x 
dt  dmn  rs  y 


n) 

n) 


transpose (new_tau_x ( : , 
transpose (new_tau_y ( : , 
reshape (transpose_tau_ 
reshape (transpose_tau_ 


:,n)); 

:,n)); 

X  ( : ,  : ,  n)  ,  1, 

y  ( : ,  : ,  n)  ,  1, 


=  zeros ([624  1679]); 

=  zeros ([624  1679]); 

=  zeros ([624  1679]); 
=  zeros  ([624  1679]); 


1679)  ; 
1679)  ; 
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%demean, detrend  the  position  thru  time 
for  n  =  1:1679 

dmn  reshape  x(:,n)  =  reshape  tau  x (:, n) -mean (reshape  tau  x(:,n)); 
dmn_reshape_y ( : , n)  =  reshape_tau_y ( : , n) -mean (reshape_tau_y ( : , n) ) ; 
dt  dmn  rs  x(:,n)  =  detrend (dmn  reshape  x(:,n)); 
dt_dmn_rs_y ( : , n)  =  detrend (dmn_reshape_y (:, n) ) ; 
n 

end% 


%create  empty  matrix 

rs_x  =  zeros ([73  23  624]); 

rs_y  =  zeros ([73  23  624]); 

new  X  =  zeros ( [23  73  624]); 
new_y  =  zeros ( [23  73  624]); 

%unwrap  and  restore  to  original  matrix  now  that  it  is 
detrended/ demeaned 
for  n  =  1:624 

rs_x(:,:,n)  =  reshape (dt_dmn_rs_x (n, :) ,  73,23); 
rs_y(:,:,n)  =  reshape (dt_dmn_rs_y (n, :) ,  73,23); 
new_x(:,:,n)  =  transpose ( rs_x (:,:, n) ) ; 
new_y ( : , : , n)  =  transpose (rs_y ( : , : , n) ) ; 
end 

%FINISH  (DETREND,  DEMEAN  U,V  COMPONENTS) 

q, 

o 

%make  v  complex 

q, 

o 

%create  empty  matrix 

complex_tau_y  =  zeros ( [23  73  624]); 

for  n=l : 624 

complex_tau_y ( : , : , n)  =  new_y ( : , : , n) *li; 
str=num2str (n) ; 
disp ( [ ' i  ' ,  str ] ) 
end 


q, 

o 

%add  u+iv 

q, 

o 

%create  empty  matrix 

complex  vector  =  zeros ( [23  73  624]); 

for  n=l : 624 

complex_vector ( : , : , n)  =  complex_tau_y ( : , : , n)  +  new_x ( : , : , n) ; 

str=num2str (n) ; 

disp ([' vector  ',  str]) 

end 
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q, 

o 

%  RESHAPE  DATA  FOR  EOF  ANALYSIS 

q, 

o 

%create  empty  matrix 

transpose_cv  =  zeros ([73  23  624]); 
reshape_cv  =  zeros ([624  1679]); 
tr_reshape_cv  =  zeros ([1679  624]); 

%reshape 

for  n  =  1:624 

transpose_cv ( : , : , n)  =  (complex_vector (:,:,n).'); 

%2 3x7 3x62 4 ->73x2 3x62 4 

reshape_cv (n, : )  =  reshape (transpose_cv (:,:, n) ,  1,  1679); 
%7 3x2 3x62 4 ->624x1 67 9 

tr_reshape_cv ( : , n)  =  ( reshape_cv (n, : ) . ' ) ; 

%624xl679  ->1679x624 
str=num2str (n) ; 
disp ([' reshape  str]  ) 
end 


q, 

o 

%create  empty  matrix  for  variance/covariance  matrix 
H  =  zeros ([1679  1679]); 

H(:,:)  =  (  (tr_reshape_cv*tr_reshape_cv' ) 7624  ); 

q, 

o 

%  FIG  Eigenvalues  and  eigenvectors. 

%  E  =  EIG(X)  is  a  vector  containing  the  eigenvalues  of  a  square 
%  matrix  X. 

%  [V,D]  =  EIG(X)  produces  a  diagonal  matrix  D  of  eigenvalues  and  a 

%  full  matrix  V  whose  columns  are  the  corresponding  eigenvectors  so 

%  that  X*V  =  V*D. 


q, 

o 

%create  the  empty  matrix  for  V  and  D 
V  =  zeros ([1679  1679]); 

D  =  zeros ([1679  1679]); 

[V,D]  =  eig(H) ; 

q, 

o 

%create  the  empty  matrix  for  the  coefficients  (PCs) 

C=zeros (1679, 624) ; 

%get  the  coefficients  of  the  eigenvector  for  each  pc  of  month  m 
%  (so  col  1  contains  all  the  coefficients  of  eigenvector  1  (1679) etc.) 

for  m=l : 624 
for  n=l:1679 

C (n, m) =V ( : , n) ' *tr_reshape_cv ( : , m) ; 
end  %n 
end  %m 

q, 

o 

%  PC(k)  time  series  k=l:1679 

%  PCI  =  Cl  Complex  in  the  form  a+ib 

C1=C(1, :) ; 

Cl  l=real(Cl);  %primary  component 
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Cl  2=imag (Cl ) ;  %secondary  component 

%  PC2  =  C2  Complex  in  the  form  a+ib 
C2=C (2, : ) ; 

C2  l=real(Cl);  %primary  component 

C2  2=imag(Cl);  %secondary  component 

%Etc . 

q, 

o 

%PLOT  PCI  time  series  (PCI  for  624  months) 
plot(l:624,  Cl_l, 

q, 

o 

%calculate  phase  time  series 
%create  empty  matrix 

phase_ts  =  zeros ([624  1]);  %1679xl679 

for  n=l : 624 

phase_ts (n)  =  [atan2 (  (imag (C(l,n))), (real(C(l,n)))  )]; 

deg_phase_ts (n)  =  rad2deg (phase_ts (n) ) ; 

end 

q, 

o 

%calculate  weighting  factor  time  series 
%create  empty  matrix 

weight_ts  =  zeros ([624  1]);  %1679xl679 

for  n=l : 624 

weight_ts (n) = [sqrt ( (real(C(l,n) ) *real(C(l,n) ) )+ (imag (C (1, n) ) *imag (C (1, n 
)  )  )  )  ]  ; 
end 

q, 

o 

%  create  the  1st  eigenvector  field 
V(:,1679) ; 
uVl=real (V ( : ,  1 ) )  ; 
vVl=imag (V ( :  ,  1 ) )  ; 

tr_uVl  =  uVl ' ; 

rs_uVl  =  reshape (tr_uVl, 73, 23) ; 
tr_rs  uVl  =  rs  uVl ' ; 

tr_vVl  =  vVl ' ; 

rs  vVl  =  reshape (tr  vVl,73,23); 
tr  rs  vVl  =  rs  vVl '  ; 

q,  _  _  _  _ 

o 

%PLOT  Eigenvector 

q, 

o 

hold  on 

quiver ( lonValues ( 1 : 2 : 73 ) ,  latValues,  tr  rs  uVl ( : , 1 : 2 : 73 ) , 

tr  rs  vVl ( : , 1 : 2 : 73 ) ) 
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o\°  o\°  o\°  o\° 


variance  calculation 

var=zeros (1679, 1) ; 
var=real (diag (D) ) 

pclvar  =  (var (1) / (trace (H) )) *100 
pc2var  =  (var (2) / (trace (H) )) *100 
pc3var  =  (var (3) / (trace (H) )) *100 
pc4var  =  (var (4) / (trace (H) )) *100 
pcSvar  =  (var (5) / (trace (H) )) *100 
pc6var  =  (var (6) / (trace (H) )) *100 

etc . 
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